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1.  Introduction 

Investigations  of  nonrigid-Earth  rotation  usually  rely  on  the  convolution  of  a 
transfer  function  with  a  rigid- Earth  solution.  A  possible  alternative  is  to  develop 
a  unified  theory  to  encompass  in  only  one  body  both  rigid  and  nonrigid  effects. 
It  is  represented  in  the  Hamiltonian  theory  by  Getino  and  Ferrandiz  (GF).  This 
theory  combines  the  convenience  of  being  a  really  analytical  theory  with  good 
practical  performance,  since  the  residuals  with  respect  to  observations  are  also 
at  the  same  ±1  mas  level  as  the  present  IERS  series. 

We  report  on  more  recent  results  we  obtained  in  the  last  year,  not  collected 
in  the  “Considerations”  paper  of  the  WGN  (joint  IAU/IUGG  Working  Group 
on  Nutations). 


2.  General  features  of  GF  theory 

The  more  relevant  features  of  the  GF  theory  are:  1)  It  is  a  unified  and  self- 
consistent  theory:  rigid  and  nonrigid  effects  are  derived  from  a  unique  Hamilto¬ 
nian  function.  2)  It  is  enterely  analytical:  Lunar  and  planetary  solutions  being 
taken  from  analytical  ephemerides  (ELP2000,  VSOP87).  The  solutions  are  ob¬ 
tained  through  a  perturbation  method  (Hori’s).  3)  Parameters  appearing  in  the 
analytical  solution  are  global,  and  can  be  determined  by  numerical  integration 
over  the  Earth’s  interior  or  by  fitting  the  series  to  observations.  4)  A  general 
Hamiltonian  formalism  is  used.  We  point  out  the  computational  advantage, 
especially  to  deal  with  second  order  (nonlinear)  perturbations.  5)  Precession 
and  nutation  can  be  treated  jointly.  Formulas  for  the  precession  of  the  nonrigid 
Earth  can  be  derived  from  the  secular  part  of  the  same  Hamiltonian. 


3.  Present  stage  of  development 

The  main  part  of  the  theory  considers  a  three-layer  Earth  composed  of  mantle, 
Fluid  Outer  Core  (FOC)  and  Solid  Inner  Core  (SIC).  Dissipation  at  the  two 
boundaries  is  modelled  similar  to  Sasao  et  al.  (1980),  taking  into  account  both 
viscous  and  electromagnetic  coupling.  Elasticity  is  treated  along  the  lines  of 
Takeuchi,  followed  by  Jeffreys  and  Vicente,  and  complemented  with  an  elastic 
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Figure  1.  Absolute  value  of  deviations  with  respect  to  IERS96.  In- 
phase  terms. 


delay  in  the  response.  Nonrigid  second  order  perturbations  (in  the  sense  of  per¬ 
turbation  theory,  that  are  quadratic  in  the  dynamical  ellipticity )  have  been  par¬ 
tially  derived.  These  terms  cannot  be  found  by  the  use  of  the  transfer  function 
approach  that  is  essentially  linear  with  respect  to  that  perturbation  parameter. 
High-frequency  terms  have  been  derived  for  the  nonrigid  model,  corresponding 
to  triaxiality  (EGS99)  and  third  and  forth  order  harmonics.  Analytical  expres¬ 
sions  for  secular  mixed  terms  are  also  already  available  (EGS2000).  Nonrigid 
contributions  from  planetary  origins  have  been  computed  in  part,  and  are  not 
negligible.  The  same  is  valid  for  the  additional  tidal  potential.  Analytical  and 
numerical  partial  derivatives  with  respect  to  all  the  parameters  are  also  available. 


4.  Accuracy  performance 

In  Figures  1  and  *2,  prograde  and  retrograde  amplitudes  of  the  main  terms  pro¬ 
vided  by  the  series  GF99-3  by  Getino  and  Ferrandiz  are  compared  with  those 
of  Mathews,  Buffet  and  Herring  1999  (MBH99),  Dehant  and  Defraigne  1997 
(DD97)  and  Scliastok  1997  (SGI  197),  showing  the  absolute  value  of  respective 
deviations  with  respect  to  IERS96.  Model  GF99-  3i  S  t  includes  the  oceanic  cor¬ 
rections  given  by  Scliastok  (1997),  while  G 199-3 ( 1)1))  uses  the  oceanic  correc¬ 
tions  provided  by  Dehant  and  Defraigne  (1997). 

Time  domain  deviations  of  GF99-3  series  with  respect  to  IER.S96  are  pre¬ 
sented  in  Figure  3  for  longitude  ( A^-sine)  and  Figure  4  for  obliquity  (Ae).  Next 
we  compare  the  time  domain  deviations  of  GF99-3  and  IERS96  with  respect  to 
the  observations  USNO9903.  Deviations  corresponding  to  our  series  GF99-3  are 
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Figure  2.  Absolute  value  of  deviations  with  respect  to  IERS96.  Out- 
of-pliase  terms. 


marked  with  a  symbol  and  those  of  the  IERS  are  marked  with  an  “o”. 

Figure  5  shows  longitude  (Ayi)  and  Figure  6  obliquity  (As). 


5.  Examples  of  new  nonrigid  contributions 

For  the  sake  of  brevity  we  present  only  some  few  terms. 

1. -  FCN  resonance  effect  on  Venus’s  direct  perturbation:  Main  terms  in 
longitude  (fias):  A tp  =  —113  cos(2 Lye  —  4 Le  +  2 h)  —  82  ain(2Xye  —  4Zjc;  +  2 h)  + 
75  sin(Xye  —  Le  +  h). 

2. -  Second  order  perturbation:  main  effect  of  FCN  resonance  on  Oppolzer 
terms  (f las ):  A tp  =  23 sin( il )  +  27 >in( 20 )• 

3. -  Tidal  potential:  main  contributions  (/. tas )  A e  =  163c.os(fl). 

4. -  Semidiurnal  nonrigid  main  contributions  (/. ias )  At/’  =  3  sin!  2o  —  2  /■  — 

2  ft)  +  4  sin(2<^). 

6.  Conclusions  and  suggestions 

Nowadays  we  are  not  forced  to  decide  on  choosing  just  one  rigid  Earth  theory  and 
one  transfer  function,  because  there  is  an  alternative  option  available,  namely 
to  choose  the  GF  unified  theory  to  encompass  in  only  one  body  both  rigid  and 
nonrigid  effects.  In  our  opinion,  future  IATJ  Resolutions  concerning  nutations 
should  be  aware  of  the  existence  of  the  unified  Hamiltonian  theory  of  the  nonrigid 
Earth,  and  not  restrict  themselves  to  the  several  rigid  Earth  theory  +  transfer 
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Figure  3.  Time  domain  deviations  of  GF99-3:  At/’  •  sine. 
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Figure  4.  Time  domain  deviations  of  GF99-3:  As. 


Figure  5.  Time  domain  deviations  witli  respect  to  observations:  Lon¬ 
gitude  (A^>). 
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Figure  6.  Time  domain  deviations  with  respect  to  observations: 
Obliquity  (Ae|, 


function  approaches.  The  reason  is  that  the  GF  theory  shows  more  capabilities 
than  any  other  method,  including  state-of-the-art  level  of  accuracy,  complete 
consistency  and  ability  to  derive  second  order  and  other  non  rigid  corrections 
near  or  up  to  the  microarcsecond  level  for  each  frequency.  The  latter  is  especially 
important  for  future  improvements,  since  they  are  not  likely  to  assume  that  the 
large  discrepancy  between  all  theories  and  observations  (measured  in  terms  of 
sigmas)  is  mainly  due  to  chaotic  dynamics  and  that  we  have  reached  the  limits 
of  predictability. 

We  have  asked  also  the  WGN  to  consider  the  possibility  of  our  Unified 
Theory  to  play  a  role  in  the  eventual  adoption  by  the  WGN  of  a  new  nonrigid 
Earth  theory.  The  requirements  of  a  new  theory  closer  to  observations  than 
IAU  4980  may  be  covered  by  the  combination  of  1)  a  main  analytical  part, 
the  Hamiltonian  theory  for  the  solid,  nonrigid  Earth  described  above,  including 
partial  derivatives  and  nonlinear  liigh-order  perturbations,  2)  high-frequency 
and  planetary  terms,  that  should  become  more  relevant  in  the  future  and  should 
be  derived  from  a  nonrigid  model  at  the  1-mas  truncation  level,  and  3)  some 
additional  corrections  corresponding,  for  instance,  to  the  more  well-known  effect 
of  outer  geophysical  fluids,  at  least  the  oceans. 

Note  this  possibility  is  not  in  opposition  to  using  successive  IERS  models  to 
fit  better  the  observations,  like  IERS’96  or  MHB2000.  Coexistence  of  a  theory 
and  a  computational  model  has  been  successful  and  stimulated  advances  for 
some  years.  In  our  opinion,  a  ‘theory’  should  provide  a  solution  as  accurate 
and  descriptive  as  possible  to  a  physical  model  and  as  close  as  possible  to  the 
actual  problem,  whereas  a  ‘computational  model’  should  provide  predictions  as 
close  as  possible  to  observations  even  though  the  underlying  physical  situation  is 
not  completely  understood  (as  happens  in  polar  motion  prediction  by  nonlinear 
dynamical  techniques ) 
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